arXiv:1506.00013vl [astro-ph.HE] 29 May 2015 


Updated Search for Spectral Lines from Galactic Dark Matter Interactions with Pass 

8 Data from the Fermi Large Area Telescope 

M. Ackermann/ M. Ajello,^ A. Albert,B. Anderson,"*^ W. B. Atwood,® L. Baldini,®’® G. Barbiellini,^’® 

D. Bastieri,®’R. Bellazzini,^^ E. Bissaldi,^^ R. D. Blandford,® E. D. Bloom,® R. Bonino,^®’^^ E. Bottacini,® 

T. J. Brandt,®® J. Bregeon,®® R Bruel,®®' R. Buehler,® S. Buson,®’®° G. A. Galiandro,®’®® R. A. Cameron,® 

R. Caputo,®’[^ M. Caragiulo,®^ P. A. Caraveo,®® C. Cecchi,^®’^® E. Charles,® A. Chekhtman,^^ J. Chiang,® 

G. Chiaro,®® S. Ciprini,23.20,24 Claus,® J. Cohen-Tanugi,®® J. Conrad,2®'26.4.27 a. Cuoco,^®' i®' i4 

S. Cutini,2®>24.20 p_ D’Ammando,28-29 A. de Angelis,®® F. de Palma,®2>3® R. Desiante,'®’®2 S. W. Digel,® 

L. Di Venere,®® P. S. Drell,® A. Drlica-Wagner,®4 C. Favuzzi,®®’ ®2 S. J. Fegan,®®" A. Franckowiak,® 

Y. Fukazawa,®® S. Funk,® P. Fusco,®®’®2 F. Gargano,®2 D. Gasparrini,2®’ 24 ,20 Giglietto,®®’®2 
F. Giordano,®®’®2 M. Giroletti,28 G. Godfrey,® G. A. Gomez-Vargas,®®’3® I. A. Grenier,®® J. E. Grove,®® 

S. Guiriec,®®’"®® M. Gustafsson,^^ J.W. Hewitt,42-43 a. B. Hill, 44 > 3,45 Horan,®®" G. J6hannesson,4® 

R. P. Johnson,® M. Kuss,®® S. Larsson,2®’26.47 ^ Latronico,®® J. Li,4® L. Li,4 ®’26 p, Longo,®"’® F. Loparco,®®’®2 
M. N. Lovellette,®® P. Lubrano,2®’ 21 D. Malyshev,® M. Mayer,® M. N. Mazziotta,®2 J. E. McEnery,®®’®® 

P. F. Michelson,® T. Mizuno,®® A. A. Moiseev,^®’®® M. E. Monzani,® A. Morselli,®® S. Murgia,®2 E. Nuss,®® 

T. Ohsugi,®® M. Orient!, 2 ® E. Orlando,® J. F. Ormes,®® D. Paneque,®^’® M. Pesce-Rollins,®® F. Piron,®® 

G. Pivato,®® S. Rain6,®®’®2 R. Rando,®’®® M. Razzano,®®’®® A. Reimer,®®’® T. Reposeur,®^ S. Ritz,® 

M. Sanchez-Conde,2®’ 2 ® A. Schulz,® C. Sgro,®® E. J. Siskind,®® F. Spada,®® G. Spandre,®® P. Spinelli,®®’®2 
H. Tajima,®®’® H. Takahashi,®® J. B. Thayer,® L. Tibaldo,® D. F. Torres,"®®’®o G. Tosti,®®’®! E. Troja,®®’®® 

G. Vianello,® M. Werner,®® B. L. Winer,®® K. S. Wood,®® M. Wood,® G. Zaharijas,®®’®® and S. Zimmer2®’2® 

^Deutsches Elektronen Synchrotron DESY, D-15738 Zeuthen, Germany 
^Department of Physics and Astronomy, Clemson University, 

Kinard Lab of Physics, Clemson, SC 29634-0978, USA 
® W. W. Hansen Experimental Physics Laboratory, 

Kavli Institute for Particle Astrophysics and Cosmology, 

Department of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 94305, USA 

^Royal Swedish Academy of Sciences Research Fellow, 
funded by a grant from the K. A. Wallenberg Foundation 
^Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and Astrophysics, 
University of California at Santa Cruz, Santa Cruz, CA 95064, USA 
^ Universitd di Pisa and Istituto Nazionale di Fisica Nucleare, Sezione di Pisa 1-56127 Pisa, Italy 
^ Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 1-34127 Trieste, Italy 
^ Dipartimento di Fisica, Universitd di Trieste, 1-34127 Trieste, Italy 
^Istituto Nazionale di Fisica Nucleare, Sezione di Padova, 1-35131 Padova, Italy 
Dipartimento di Fisica e Astronomia “G. Galilei”, Universitd di Padova, 1-35131 Padova, Italy 
^^Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 1-56127 Pisa, Italy 
Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 70126 Bari, Italy 
Istituto Nazionale di Fisica Nucleare, Sezione di Torino, 1-10125 Torino, Italy 
Dipartimento di Fisica Generale “Amadeo Avogadro” , 

Universitd degli Studi di Torino, 1-10125 Torino, Italy 
^^NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 
Laboratoire Univers et Particules de Montpellier, 

Universite Montpellier, CNRS/IN2P3, Montpellier, France 
Laboratoire Leprince-Ringuet, Ecole poly technique, CNRS/IN2P3, Palaiseau, France 
^^Consorzio Interuniversitario per la Fisica Spaziale (CIFS), 1-10133 Torino, Italy 
INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, 1-20133 Milano, Italy 
Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 1-06123 Perugia, Italy 
Dipartimento di Fisica, Universitd degli Studi di Perugia, 1-06123 Perugia, Italy 
22 College of Scienee, George Mason University, Fairfax, VA 22030, 
resident at Naval Research Laboratory, Washington, DC 20375, USA 
^^Agenzia Spaziale Italiana (ASI) Seience Data Center, 1-00133 Roma, Italy 
^"^INAF Osservatorio Astronomico di Roma, I-OOO 4 O Monte Porzio Catone (Roma), Italy 
Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden 
The Oskar Klein Centre for Cosmoparticle Physics, AlbaNova, SE-106 91 Stockholm, Sweden 
The Royal Swedish Academy of Sciences, Box 50005, SE-IO 4 05 Stockholm, Sweden 
^^INAF Istituto di Radioastronomia, 40129 Bologna, Italy 
Dipartimento di Astronomia, Universitd di Bologna, 1-40127 Bologna, Italy 


2 


Dipartimento di Fisica, Universitd di Udine and Istituto Nazionale di Fisica Nucleare, 

Sezione di Trieste, Gruppo Collegato di Udine, 1-33100 Udine 

Universitd Telematica Pegaso, Piazza Trieste e Trento, 48, 80132 Napoli, Italy 
Universitd di Udine, 1-33100 Udine, Italy 

Dipartimento di Fisica “M. Merlin” delTUniversitd e del Politecnico di Bari, 1-70126 Bari, Italy 
Center for Particle Astrophysics, Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 
Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan 
^^Istituto Nazionale di Fisica Nucleare, Sezione di Roma “Tor Vergata”, 1-00133 Roma, Italy 

Departamento de Fisica, Pontificia Universidad Catolica de Chile, 

Avenida Vicuna Mackenna 4860, Santiago, Chile 

Laboratoire AIM, CEA-IRFU/CNRS/Universite Paris Diderot, 

Service d’Astrophysique, CEA Saclay, 91191 Gif sur Yvette, France 

Space Science Division, Naval Research Laboratory, Washington, DC 20375-5352, USA 

'^^NASA Postdoctoral Program Fellow, USA 

Georg-August University Gottingen, Institute for theoretical Physics - Faculty of Physics, 

Friedrich-Hund-Platz 1, D-37077 Gottingen, Germany 

Department of Physics and Center for Space Sciences and Technology, 

University of Maryland Baltimore County, Baltimore, MD 21250, USA 

Center for Research and Exploration in Space Science and Technology 
(CRESST) and NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA 

School of Physics and Astronomy, University of Southampton, Highfield, Southampton, S017 IBJ, UK 

Funded by a Marie Curie lOF, FP7/2007-2013 - Grant agreement no. 275861 

Science Institute, University of Iceland, IS-107 Reykjavik, Iceland 

Department of Astronomy, Stockholm University, SE-106 91 Stockholm, Sweden 

Institute of Space Sciences (lEEC-CSIC), Campus UAB, E-08193 Barcelona, Spain 

Department of Physics, KTH Royal Institute of Technology, AlbaNova, SE-106 91 Stockholm, Sweden 

Department of Physics and Department of Astronomy, 

University of Maryland, College Park, MD 20742, USA 

Hiroshima Astrophysical Science Center, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan 

Center for Cosmology, Physics and Astronomy Department, 

University of California, Irvine, CA 92697-2575, USA 

Department of Physics and Astronomy, University of Denver, Denver, CO 80208, USA 

Max-Planck-Institut fiir Physik, D-80805 Munchen, Germany 

Funded by contract FIRB-2012-RBFR12PM1F from the Italian Ministry of Education, University and Research (MIUR) 

^^Institut fiir Astro- und Teilchenphysik and Institut fiir Theoretische Physik, 
Leopold-Franzens-Universitdt Innsbruck, A-6020 Innsbruck, Austria 

Centre d’Etudes Nucleaires de Bordeaux Gradignan, IN2P3/CNRS, 

Universite Bordeaux 1, BP120, F-33175 Gradignan Cedex, France 

^^NYCB Real-Time Computing Inc., Lattingtown, NY 11560-1025, USA 
Solar-Terrestrial Environment Laboratory, Nagoya University, Nagoya 464-8601, Japan 
Institucio Catalana de Recerca i Estudis Avangats (ICREA), Barcelona, Spain 

Department of Physics, Center for Cosmology and Astro-Particle Physics, 

The Ohio State University, Columbus, OH 43210, USA 

Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 
and Universitd di Trieste, 1-34127 Trieste, Italy 

Laboratory for Astroparticle Physics, University of Nova Gorica, Vipavska 13, SI-5000 Nova Gorica, Slovenia 

(Dated: June 2, 2015) 



3 


Dark matter in the Milky Way may annihilate directly into 7 rays, producing a monoenergetic 
spectral line. Therefore, detecting such a signature would be strong evidence for dark matter 
annihilation or decay. We search for spectral lines in the Fermi Large Area Telescope observations 
of the Milky Way halo in the energy range 200 MeV - 500 GeV using analysis methods from our 
most recent line searches. The main improvements relative to previous works are our use of 5.8 
years of data reprocessed with the Pass 8 event-level analysis and the additional data resulting 
from the modified observing strategy designed to increase exposure of the Galactic center region. 
We searched in five sky regions selected to optimize sensitivity to different theoretically-motivated 
dark matter scenarios and find no significant detections. In addition to presenting the results from 
our search for lines, we also investigate the previously reported tentative detection of a line at 133 
GeV using the new Pass 8 data. 


I. INTRODUCTION 


Cosmological observations reveal that ~80% of the matter in the Universe is dark matter (DM). Cosmic microwave 
background measurements, galactic rotation curves, gravitational lensing (among others) each provide strong evidence 
for the existence of DM m- One of the leading DM candidates are the Weakly Interacting Massive Particles (WIMP), 
though other candidates such as gravitinos may also account for some, or all of the observed DM (see Ref. mis] for 
recent reviews of particle DM candidates). 

WIMPs and other DM candidates may annihilate or decay to Standard Model particles, which would produce 7 
rays. Gamma-ray signals from DM annihilation or decay are expected to typically produce a broad spectral signature 
(see Ref. IS [7] for recent indirect DM search reviews). The difficulty in detecting such a signal lies in distinguishing 
it from other standard astrophysical processes. However, if the DM annihilates or decays into a photon and a neutral 
particle (such as another photon or Z boson), approximately monoenergetic 7 rays will be produced in the rest frame. 
For non-relativistic DM particles, this would give rise to a monoenergetic photon signal in the otherwise smooth 
spectrum of the standard astrophysical emission. Such a sharp spectral signature is not expected from standard 
astrophysical mechanisms, though non-DM induced mechanisms have been proposed [e.g.,| 8 j. The branching fraction 
of monoenergetic DM annihilation channels is typically loop-suppressed, ~ 10“^(cru) — 10“^(cm), where (tru) 

is the total velocity-averaged DM annihilation cross section and {o'v)^^ is the cross section for DM annihilation to 
two 7 rays iHn]. The total annihilation cross section that would produce the currently observed abundance of DM 
in the Universe is {av) « 3 x 10“^® cm“^s“^, assuming that DM is a thermal relic |I4) . 

In this paper, we search for 7 -ray spectral lines using data obtained by the Large Area Telescope (LAT) [T^ on board 
the Fermi 7 -ray Space Telescope. The LAT has been surveying the 7 -ray sky in an energy range of 20 MeV to over 
300 GeV since 2008, making it an ideal instrument to search for 7 rays from DM interactions. WIMP DM candidates 
typically have ^GeV to ~TeV scale masses; therefore, WIMP annihilations would produce 7 rays detectable by the 
LAT. Specifically, DM annihilations directly into pairs of 7 rays will create a spectral line at the rest mass energy of 
the DM particle. Besides being a strong indication of WIMP interactions, spectral lines are one of the best ways to 
search for gravitino (another DM candidate) decay [16] . 

The LAT Collaboration has published four previous line searches, each improving and expanding the analysis relative 
to the previous works. Our first analysis searched for lines from 30 to 200 GeV using 11 months of data [T7]. The 
search using two years of data expanded the energy range to search for lines down to 7 GeV and made use of control 
regions to estimate systematic uncertainties [H]. The 3.7 year analysis covered the energy range from 5 to 300 GeV 
and included a thorough investigation and quantification of the systematic uncertainties |19j . That work also included 
a detailed investigation of a tentative line-like feature near 133 GeV reported in the region around the Galactic center 
(GC) [201 mi- An analysis extending to lower energies (100 MeV to 10 GeV), where systematic uncertainties typically 
dominate, was performed using 5.2 years of data [22] . In that paper, the systematic uncertainties were quantified, 
and also incorporated in the likelihood analysis. These past analyses used previous iterations (or ‘Passes’) of the LAT 
event reconstruction and classification [II1I23I. 

In this work, we update the results from our previous line searches using 5.8 years of Pass 8 data. We present an 
analysis performed over three decades in energy that incorporates systematic uncertainties estimated using control 
regions. Section [nl Sec.|IIH and Sec. |IV| outline the event selection, regions of interest, and energy dispersion modeling 
respectively. Section|V|describes the methods for fitting the 7 -ray spectrum and Sec.|V^details the systematic studies. 
Section |VII| presents the results of searching for a monoenergetic line in the 7 -ray spectrum. Section VIII goes into 
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greater detail on the previously reported tentative detection near 133 GeV in the region around the GC. Finally, 
Sec. IIXI summarizes the results. 


II. DATA AND EVENT SELECTION 


The LAT Gollaboration has recently developed an extensive update of the event reconstruction and classification 
(Pass 8) which improves the performance of the LAT and reduces systematic uncertainties |24j. Pass 8 events are 
classified according to the estimated accuracy of the direction reconstruction. Each event was assigned a PSF ‘type’ 
from PSFO (worst) to PSF3 (best), with each type containing about one quarter of the total number of events in 
logarithmic energy bins. Events are also classified according to the estimated accuracy of the energy reconstruction. 
Each event is assigned an energy dispersion ‘type’ (ct) from EDISPO (worst) to EDISP3 (best) according to the 
reconstruction quality estimator ‘BestEnergyProb’. Each energy dispersion event type contains roughly one quarter 
of the events in log(E/MeV) energy bins. EigureH (left) shows the energy dispersion DeS ((A' — E)/E) at 100 GeV 
for each EDISP type. The energy resolution ( 68 % containment) is shown in Fig. (right). The expected energy 
resolution is greatly improved in EDISP3 compared to EDISPO. 




FIG. 1. P8_CLEAN_V5 energy dispersion at 100 GeV (left) and energy resolution as a function of energy (right) for the four 
event types. E is the true energy of the generated 7 ray and E' is the reconstrncted energy. 

Figure 1^ compares the energy dispersion in Pass 8 and Pass 7REP for the energy range 5 < log(E)/MeV < 5.2. 
While the energy resolution has not significantly improved in Pass 8 , the number of events accepted in each event 
class is higher. The increase in acceptance is due, in part, to the ability of event reconstruction in Pass 8 to isolate 
the gamma-ray events when a cosmic ray is in near coincidence with the gamma ray, more generally referred to as 
pile-up or ‘ghost’ events [53]. The largest impact relevant to this analysis is an increased effective area in Pass 8 by 
^ 30% for events above 10 GeV for the Glean class. We note that the specific ‘Clean’ event selections are different 
between Pass 7REP and Pass 8 . However, both are defined such that the resulting residual CR rate is at or below 
the extragalactic 7 -ray background above 100 MeV. 

We used data from a 5.8 year period with the P8_CLEAN event selection to search for spectral lines in the energy 
range 200 MeV to 500 GeV. To allow for sidebands in each fit, we used LAT 7 -ray data from 100 MeV to 750 GeV. 
Details on additional event selections are given in Tab. |l| We split the data into Celestial and Earth Limb datasets 
using cuts on both the instrument rocking angle (6^) and the event zenith angle (Oz); see Fig. 1 of Ref. [T3| for a 
schematic depicting 6r, Oz-, and event incident angle 9. Gamma rays from the Earth’s Limb are produced by cosmic-ray 
(CR) interactions in the Earth’s upper atmosphere so it is used as a control region where no true spectral lines are 
expected from dark matter annihilation or decay (see App. 0 . 

We extract our signal regions of interest (ROIs) and other control regions from the Celestial data set, where we 
require 9z < 100° in order to remove emission from the bright Earth’s Limb. No point source masking is applied since 
at lower energies the width of the point spread function (PSF) increases significantly. Using an energy-dependent 
point source mask, as was done in our 3.7-year analysis, would result in removing most of the dataset below ~1 GeV. 
The initial data reduction and all of the exposure calculations were performed using the LAT Science Tool.^ version 
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The ScienceTools and documentation are available at http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.htinl 
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FIG. 2. Normalized energy dispersion in the range 5 < log(E)/MeV < 5.2 from MC with CLEAN event class selection for Pass 
8 and Pass 7REP respectively. E is the true energy of the generated 7 ray and E' is the reconstructed energy. 


09-33-03, and the P8_CLEAN_V5 instrument response function^ 


TABLE I. Summary table of data selections. 


Selection 

Celestial 

Earth Limb 

Observation period 

August 4, 2008 - April 30, 2014 August 4, 2008 - April 30, 2014 

Mission Elapsed Tim^ (s) 

[239557417-420595073] 

[239557417-420595073] 

Energy Range (GeV) 

[0.1-750] 

]2.5-750] 

Rocking Angle {6r) 

< 52° 

> 52° 

Zenith Angle {dz) 

< 100 ° 

111° kOz < 113° 

Data qualitj|^ 

Yes 

Yes 


Mission Elapsed Time is the number of seconds since 00:00:00 UTC January 1, 2000. 
DATA.QUAL == 1 && LAT.CONFIG ==1 


III. REGIONS OF INTEREST 

Within our Celestial dataset, we define our signal ROIs to be the same as those considered in Ref. [TH]. These 
ROIs were optimized for either annihilating or decaying DM and for different profiles of the spacial distribution of 
the DM distribution in the Galaxy. We used four models of the DM distribution: Navarro-Frenk-White (NFW) |25) . 
an adiabatically contracted NFW (NFWc), Einasto [26], and a cored isothermal profile [27] . 

A generalized NFW profile [IS] is given by 


p(r) = - — - ( 1 ) 

^ (r-/r,)7(l + r/r,)3-7 

with Ts = 20 kpc. The NFW and NFWc correspond, respectively, to the cases where 7 = 1 and 7 = 1.3. We also use 
an Einasto profile defined as: 


p{r) = poexp{-( 2 /Q;)[(r/r«)“ - 1 ]}, 


( 2 ) 


^ The IRF used (P8_CLEAN_V5) is not one that will be released publicly. P8_CLEAN_V6 will be the first released. Both have been derived 
with the same data set and have only minor technical differences that do not affect this analysis. 
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where rg = 20 kpc and a = 0.17 [26j . Finally, we consider a cored, isothermal profile given by: 


p{r) 


Po 

1 + (r/rg)2 


(3) 


with Tg = 5 kpc. We normalize all profiles by fixing the local DM density to pijQ = 8.5 kpc) = 0.4 GeV cm ^ 
Figure [^compares these DM distributions. 



FIG. 3. The four DM profiles as a function of the distance from the GC, r. 

Each ROI is defined as a circular region centered on the GC. We mask the Galactic plane (GP) except for a 12° x 10° 
box centered on the GC. The ROIs are named after the value of i?GC they subtend (e.g., for R3, Roc = 3°). The 
annihilation ROIs are R3, R16, R41, and R90 optimized for the NFWc, Einasto, NEW, and Isothermal DM profiles 
respectively. R180 is the optimal ROI to search for spectral lines from DM decay (e.g. % —?► v^). In the GP, longitudes 
further than 6° from the GC are removed from all ROIs larger than R3. This is because we do not expect a large 
DM signal in that region and the y-ray emission is dominated by standard astrophysical sources. For details on the 
ROI-optimization procedure, see App. B of ref m- 

Figure]^ shows the counts map of the Celestial dataset in R180 with the outlines of the other signal ROIs, and the 
exclusion of the GP. We note that these ROIs are different from those in Ref. [32] , which optimized the ROIs only at 
lower energies where the fits are dominated by systematic uncertainties. 

For the smallest ROI, R3, the effects of leakage both in and out of the ROI from the PSF dominate the search 
region. The 68% containment radius of PSF3 is ^ 2° at 200 MeV and improves to ~ 0.4° at I GeV. Therefore, we 
only use PSF3 events for fits in R3 below 1 GeV. For all other fits, events from every PSF type are used. 


IV. ENERGY DISPERSION MODELING 

The energy dispersion is the probability density of reconstructing a true energy (E) as E'. The LAT Collaboration 
parametrizes the energy dispersion probability density function of E' in cos6 and true energy. In this work, we define 
an effective energy dispersion {D^s), where we have averaged over cos9 assuming an isotropic sourc^ 

As in Ref. [12], we find D^s for a given E is well-described as a sum of three Gaussians. We create separate 
Deff models for each of the EDISP types. Specifically, this is done as follows: 


Defi{E'; E, et) 


E 




(4) 


^ We note that values ranging from 0.2 - 0.85 GeV cm“^ are possible at the present [29t431| . Assuming a different local DM density would 
simply scale our limits for DM annihilation and decay by a factor inversely proportional to p(rQ)^ or ^(^0) respectively 
^ Our effective energy dispersion Deff can be evaluated using the publicly available LAT energy dispersion models by performing an 
acceptance-weighted averaging along cos9. 
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FIG. 4. Counts map for the Celestial dataset binned in 1° x 1° spatial bins in the R180 ROI, and plotted in Galactic 
coordinates using the Hammer-Aitoff projection. The energy range is 1-750 GeV. Also shown are the outlines of the other 
ROIs (R3, R16, R41, and R90) used in this search. The GP region with longitude greater than 6° from the GC and latitude 
smaller than 5° is removed from all signal ROIs. 


where 03 = 1 — 02 — oi (with Ui > 0 required) and cti > CT 2 > ( 73 . We fit the triple Gaussian model at energies from 
100 MeV < E < 1 TeV in logarithmic steps of 0.25. Then we can define I?eff for any energy by interpolating the 
parameters of the Gaussian. This method differs slightly from that in Ref [19] by using EDISP type as the second 
variable (in addition to E) in the “2D” model instead of P^- By modeling the energy dispersion separately for 
each EDISP type, we are able to give higher weight to events with a better energy reconstruction. Using the EDISP 
types adds extra information in the fit and improves the statistical power over a “ID” model by ~ 10-15% depending 
on energy. 


V. FITTING 
A. Fitting Procedure 

To fit for spectral lines, we use a maximum likelihood procedure in sliding energy windows in each of the five ROIs 
described in Sec. |III[ We fit at a fixed E^ at the center of the energy window. We increment E^ in steps of 0.5 
asiEj), where aslE^) is the energy resolution (68% containment) of the LAT at E^^. We perform our fits in the 
energy domain and define both a background spectrum model (Gbkg) and a signal spectrum model (Gsig). We do 
not incorporate spatial information in our fits since it would make the resulting flux limits dependent on the DM 
distribution profile assumed. Rather, we perform a generic search for monoenergetic signals in each ROI. Since we 
fit in narrow energy windows, we approximate the gamma-ray background from diffuse and point sources as a simple 
power law. The resulting expected distribution of counts is: 


Gbkg(i?'|rbkg,nbkg) = a S{E'), (5) 

where Tbkg is the power-law index, Eq is a reference energy set to 100 MeV, and £{E') is the energy-dependent 
exposure averaged over each ROI, which is needed since the fit is performed in count space. The normalization factor 

a is defined such that the total number of background events in the fit window is ribkg = f cn £{E')dE'. 

We did not explicitly convolve our background model with the energy dispersion (i.e. for Cbkg we assume E'= E). 
For fits below 200 MeV, this approximation is not valid and significantly degrades the goodness of the fits. Therefore, 
we limit our search range to E.^ > 200 MeV. 

Our signal spectrum is Csig{E'\E^) = n'^^^Defi(E'\E^). We account for systematic uncertainties that may induce a 
false line-like signal or mask a true line-like signal in our fitting by using the procedure described in Ref. [22j . This 















is especially important for fits with very small statistical uncertainties (see below). We include a nuisance parameter 
by treating the best-fit number of signal events (nbg) as the sum of the true number of signal events (nsig) and a 
systematic offset (n-syst) such that nsig= ?^syst- We constrain Usyst by modeling it as a Gaussian with a fixed 
width tJsyst with zero mean, where cTsyst = i^/syst x ^eff and is determined based on fits in control regions (see Sec. VII. 

Our overall model to fit for a line at is: 


where 


C{E'\0) 


(^sigDes{E'\E^) + a 



X Ggyst (rigyst , l>0ff) , 


( 6 ) 


Gsyst(^syst; ^eff) 





(7) 


Des{E'\Ery) is a weighted sum over the four EDISP type^ and 0 represents the model parameters E^, Pbkg, n'^ig, 
and ribkg- 

Unbinned fits are performed when there are fewer than 10,000 events in the fit window. When there are more 
than 10,000 events, we bin the data in 63 energy bins to avoid large computation times. Each energy bin is much 
narrower than the energy resolution of the LAT, making the binned fits a close approximation to the unbinned fits, 
which we confirmed in several test cases. We fit for a monoenergetic signal, and our results are applicable to any 
sharp spectral feature much narrower than the LAT energy resolution. We will interpret our results in the context of 
DM annihilation or decay in the next section. 

To discuss uncertainties involved in a line search, we dehne a quantity called the fractional signal (/) |19ll22j . which 
can be thought of as the fractional size of a line-like signal around the peak. Specifically / is 


/ — IT'sig/beS , ( 8 ) 

where beg is the effective background below the signal peak. The number of effective background counts in a given 
energy window [E~ ,E0] is calculated as: 


beS = 


N 


i^(rbkg)bkg,fc / 


(9) 


where the summation runs over 63 energy bins in each fit window, N is the total number of events in the fit, 
Fgig fc and Ebkg.fc are the binned probability distribution functions for the signal and background models: Egig ^ = 
i^dE' Csig/usig and Ebkg.fc = ).dE' Cbkg/?T-bkg, respectively. This definition of bgg is different than that 

presented in Sec. VIA of Ref. [H]. However both approximate the number of background events under the peal]^ 
By expressing &eff in terms of a sum of a function containing the signal and background models, it is straightforward 
to calculate this quantity for analyses beyond our line search. One could also expand the sum to include spatial bins 
as was done in a recent search for a DM signal in the Large Magellanic Cloud [33]. This quantity is especially useful 
when describing systematic uncertainties since rigig and beS for line-like features induced by systematic uncertainties 
are both expected to scale as N. The number of events in the energy window varies greatly across our energy ranges 
and ROIs. In RI80, the fit for E^= 214 MeV has 38.9 million events, while the fit at E^= 467 GeV in R3 has 52 
events. The quantity / in effect normalizes the large differences in the numbers of events. Note that the significance 
of a systematically induced line-like feature at a fixed / will scale as VW. 

Depending on the energy window, the fit will be dominated by either statistical ((5/stat = 1/V^eff) or systematic 
(dfsyst) uncertainties. It is critical to account for the systematic uncertainties in the fitting procedure since small 
line-like features (/ ~ 0.01) can be statistically significant, if Sfsyst is not accounted for, when Sfstat ^ <5/syst- For 
example, our fit at 947 MeV in R41 would have a local significance (see definition below) of 8.9 (t if Sf syst were 
neglected. However, in that case / = 0.008, which is well within the systematic uncertainty range (see Sec. |VI| ) and 
therefore cannot be deemed a detection. 


® We weight using the observed distribution of EDISP types in each fit ROI and energy window. 

® Eqn. (j^ approximates the statistical uncertainty on n^ig such that (Srigig si l/V^’efE (see Sec 3.1 of Ref. 1320 . 
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If we increase the energy window width, the statistical uncertainty decreases because more events are introduced 
into the fit. However, the systematic uncertainty (Jrisyst = <5/syst x ^eff) would then increase because the power-law 
approximation of the background energy spectrum becomes less valid. In our previous works we characterized the 
window width in terms of as (68% containment) at E^. The 3.7-year search [19) . for which most fits were dominated 
by statistical uncertainties, used ±6(T£; windows. The 5.2-year low-energy search [52] was dominated by systematic 
uncertainties and we found ±2o’£; windows to be optimal. In order to have a common definition across our entire 
energy range, we chose to define the window width to vary with E-f-. ±0.5E~f. This causes the fit range to decrease 
naturally relative to (Te{E^) at low energies and widen relative to (Je{E^) at higher energies. A common window 
definition is useful since the size of line-like features observed in our control regions (which determines 5/syst) varies 
depending on the window width. 


B. Signal Significance and Trials Factor 


We define the local significance of a fit as the square root of the test statistic (siocai = ’s/TS), where the TS is 
defined as twice the logarithm of the ratio between the likelihood maximized for the signal hypothesis and for the null 
hypothesis (rigig = 0): 


TS = 21n 


Ei^Tlgig — ^Tsig.best) 
>C(nsig = 0) 


( 10 ) 


Since we perform many fits (121 fit energies in five ROIs), the local significance (siocai) must be corrected by a 
trials factor to obtained the global significance (Sgiobai)- Since our fits largely overlap spatially and in energy, simply 
assuming all 605 fits are independent would severely overestimate our trials factor. Therefore, we determine the trials 
factor using background-only simulations. We follow the procedure outlined in Section VB of Ref. m- We created 
1000 pseudo-experiments with our full-fit energy range for our five ROIs. The pseudo-experiments had a power-law 
energy spectrum with a rbkg= 2.3. We created five independent data sets at each energy which correspond to the 
spatially independent pieces of our nested ROIs. The number of events in each piece were taken from a Poisson 
randomization of the number of events seen in the actual data. We then combined the appropriate subsets to create 
the complete simulated data set for the whole region under study. For each pseudo-experiment (i.e. a simulation of 
one full search across the entire energy range and in all ROIs), we find the largest local significance (Smax)- We can 
then use the cumulative distribution of Smax to determine Sgiobai for a given siocai- We show both the Smax distribution 
and the siocai to Sgiobai conversion for our line search in Fig. We empirically found that the effective number of 
independent trials is 353 ± 11. In each of the five ROIs, the best fit number of independent fits is ^ 72.6 taking into 
account the overlap in neighboring energy windows. 




FIG. 5. Derivation of the global significance for a given local significance in this line search. Based on 1000 background-only 
pseudo-experiments of our full fit range. (Left) Distribution of the largest siocai value (smax) obtained in each pseudo-experiment. 
Each pseudo-experiment consisted of 605 fits. (Right) Corresponding sio^ai to Sgiobai conversion. 
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VI. SYSTEMATIC UNCERTAINTIES 

To estimate the systematics, we perform scans for lines in control regions where the background events vastly 
outnumber the signal events. Our main control region is the GP excluding the GC (the white region in Fig. |^. Any 
line-like features detected along the GP would be induced by systematic uncertainties such as modeling deficiencies. 

The two largest sources of systematic uncertainties are our modeling the background flux spectrum as a power 
law and our approximation of the energy-dependent variations in the exposure. Since we do not mask point sources, 
the Galactic diffuse emission and point sources are considered together and approximated as a power law in each 
narrow energy window. Also, as was discussed in Ref |19j . fine energy-dependent variations in the LAT effective 
area are difficult to model accurately. Any discrepancy in the overall background model at (see Eqn. ([^) can be 
compensated by an excess or absorption-like feature. We assume that the level of systematic uncertainties observed 
in the GP will also be present in our signal ROIs. 

Following Ref. [52], we quantify the level of the systematic uncertainties by their fractional size d/syst (see Sec. M 
for details on /). In this way we can directly apply the Sfsyst observed in the GP to our signal ROIs. We scan Iot 
spectral lines from 200 MeV to 500 GeV in 31 10° x 10° boxes along the GP. Specifically, we scan in regions where 
|5| < 5° and I > 35° or I < 325°. The results from this scan are shown in Fig. 



E, (MeV) 


FIG. 6. Fractional signals (/) in the GP. Dots show observed / in 31 10° x 10° boxes along the GP. The solid red line is the 
average of the statistical uncertainties of the individual boxes. The blue dashed line is the value we chose to characterize <5/syst 
from modeling biases; see text for details. 


The general trend of the fractional signal with energy shown in Fig. are due to the interplay between the 
background model and energy-dependent variations in the exposure, while the spread at each energy is due to variations 
from region to region in the astrophysical emission. To further investigate the energy-dependent variations in the 
exposure, we analyzed two other control data sets: the Vela Pulsar and the Earth Limb (see App. Note that the 


general energy-dependent behavior of the observed fractional signals in each control region differs (see Fig. 14). This 


is due to differences in the deficiencies of our background models for each control region compared to the observed 
emission. 

Figure ^ shows that the fits start to become dominated by statistical uncertainties (d/stat > <5 /syst) around 
Ej = 6 GeV. Therefore, we set (5/gp = 0.015 as the systematic level from the GP since it is the 68% contain¬ 
ment of /gp for Ej < 5 GeV. We note that the transition to dominance by statistical uncertainties is dependent on 
the total number of events in the fit, not E^. For Sfsyst — 0.015 this occurs around N ^ 10,000. Additionally, the 
value of Sfsyst we find h ere is slightly higher than that used in our 5.2-year search [22| since we are using wider energy 
windows (see Sec. VA). 


An additional potential source of systematic uncertainty is cosmic-ray contamination. To estimate SfcR, we use 
the method described in Refs dang. Though a subdominant effect, cosmic-ray contamination is most significant 
at high Galactic latitudes. To estimate this uncertainty, we perform fits for lines using a sample of events in R180 
that contains a high cosmic-ray contamination — a background enriched (or “dirty”) sample. These are events that 
pass the P8_S0URCE selection, but not the P8_CLEAN selection. We then find the 68% containment of the / values 
observed in our background enriched sample to define d/cR,dirty To obtain the appropriate 6 fen value to use in our 
P8_CLEAN signal dataset, we scaled (5/cr, dirty using the y-ray acceptance ratio between P8_CLEAN and P8_SDURCE and 
the observed number of events in both event selections. See App. D5 in Ref [IS] for more details. We estimate SfcR ^ 
0.01 in R180 and R90. In R41, R16, and R3, SfcR is negligible (< 0.003). Adding SfcR in quadrature with Sfop 
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gives the total systematic uncertainty in each ROI: 6 fsyst= 0.016 for R180 and R90 and 6 fsyst= 0.015 for R41, R16, 
and R3. 


VII. FITTING RESULTS 

In our search for y-ray spectral lines in five ROIs, the data did not yield any globally significant lines. Our 
most significant fit occurred at 115 GeV in R16 and had a local significance of 2.8cr, which corresponds to a global 
significance of O.dcr (see Fig.[^. In the case of a null result, the local fit significance will follow a one-sided Gaussian 
according to Ghernoff’s theorem [31]: half the local fit significances will be zero since we require Ugig > 0. In Fig. 
we fit our local significance distribution to a one-side Gaussian function and find a best-fit width of 0.81 ± 0.03. This 
is close to the expected value of one; however, a width less than one suggests our Sfsyst determined a priori is a slight 
overestimation. 



FIG. 7. Distribution of local significance values for all our fits. Best fit to a one-sided Gaussian is shown. 


Since no significant spectral lines are detected, we set flux upper limits for monoenergetic emission. We also derive 
limits for DM annihilation and decay from our flux limits. We obtain a 95% confidence level (GL) counts upper limit 
by increasing rigig until the logarithm of the likelihood decreased by 1.36 (2.71/2) with respect to the maximum. Using 
the average exposure in each ROI at (SjioiiE^)), we can find the 95% CL monoenergetic flux upper limit using: 


If the spectral line is produced by DM annihilation directly into a pair of 7 rays the expected differential flux is given 
by: 


( 12 ) 


where todm is the mass of the DM particle, ( ) = 26{Ey — E'), and E^ = rriDM- For lines produced by DM 

V / ann 

decay into a 7 ray and a second neutral particle, the expected flux is given by: 


IIe 


decay 


I 1 1 fdN-y 

dTT tdm rnoM \ dE 


rROI 


dJ, 


decay * 


decay 

c?D 




(13) 


where is the DM lifetime, ( 7 ^) = 5{E^ — E'), and E^ = The “J factors” (Jann/decay) are 

V / decay 

proportional to the expected intensity of 7 -ray emission from DM annihilation or decay in a given ROI assuming a 
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specific DM density distribution p{r). They are defined as an integral over the line of sight of the DM density: 


dJa, 


dVL 


ds p{rY 


(14) 


and 


^Jdecay 

dU 


ds p(r). 


(15) 


f l.O.S. 


Using the 
upper limits on 
Tabs. UHlYl 
Figure ^ shows the 


,„o 95% CL upper limits derived using Eqn. 0, we solve Eqn. ( [T^ and Eqn. ( |I^ for the 95% CL 
av)^j and the 95% CL lower limits on tum respectively. The flux limits and DM limits are given in 


{av)^j 95% CL upper limits in our four ROIs optimized for sensitivity to DM annihilation 


and Fig. ^shows the t^m lower limits in R180. Also shown are the corresponding limits from our previous 3.7-year 
analysis ^ and our previous 5.2 year analysis [52]. Two main factors contribute to the differences in these three sets 
of limits: different depths of exposure, and different approaches for the treatment of systematic uncertainties. As was 
discussed in Sec. [IT] while the acceptance of the LAT increased in Pass 8, the energy resolution did not significantly 
improve. The results for each ROI benefited from the increased exposure due to the larger effective area in Pass 8. 
Also, our smallest ROIs (R3 and R16), benefited from the increased exposure of the GC region during the 6th year 
of data taking: from December 4th, 2013 to December 4th 2014, Fermi operated in a modified observing mod^that 
roughly doubled the rate of increase of exposure in the GC relative to normal survey mode. 

The 3.7-year analysis did not incorporate systematic uncertainties into calculating the limits. As was shown in 
Ref 1^, accounting for systematic uncertainties makes the results more robust, especially for fits with a large number 
of events where the systematic uncertainties dominate. In our 5.2-year analysis, we chose a conservative Sfgyst value 
that resulted in all of the fits having a local significances less than Icr. In this work, we used a more realistic (5/syst 
value, which results in a distribution of the local fit significances that is significantly closer to a one-side Gaussian 
function (see Fig. [^. Therefore, on average, our current limits should represent a greater improvement over the 
5.2-year results than would be expected solely from the increased exposure, since the 5.2-year analysis was more 
conservative. 


VIII. THE LINE-LIKE FEATURE NEAR 133 GEV 


The 7 -ray spectrum in the energy range near 133 GeV has been of particular interest after a potential signal was 
reported based on 3.7 years of data for a small ROI containing the GC region Ennu. A similar, yet not globally 
significant, feature was also reported by the LAT Collaboration m- There have been two relevant developments 
since the previous results. The first is greater exposure toward the GC, due in part to the modified observing 
strategy described in Sec. |VIT1 The second is the implementation of the Pass 8 event classification that we use here 
(see Sec. 0. In this section we first compare the overlapping events between Pass 7REP and Pass 8 (Sec. VIIIA I, 


and then measure how the apparent signal has evolved with additional data in both R3 (Sec. VIIIB) and the Earth 
Limb (Sec. VIIIC). 


A. Event-Level Comparison of Pass 7REP and Pass 8 


We first compare the reconstructed energies of events that are in both the P7REP_CLEAN and P8_CLEAN event classes 
(or event selections) and in the smallest ROI (R3). Events must pass the event selections as outlined in Tab. [ij 
be located in R3, and must have measured energies in Pass 7REP greater than 20 GeV to be considered for this 
comparison. A comparison of these events and events with energies between 120-150 GeV is shown in the left panel 
of Fig. 10 The distribution of energy differences between Pass 7REP and Pass 8 is similar in shape for events in the 


window around 133 GeV as for all events above 20 GeV. We made similar studies in the ROIs outlined in Sec. m in 
the Earth Limb, and found similar results. 


^ http://fermi.gsfc.nasa.gov/ssc/proposals/alt_obs/obs_modes.html 
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FIG. 8. 95% CL {av)ryry uppci limits for each DM profile considered in the corresponding optimized ROI. The upper left panel 
is for the NFWc ( 7 = 1 . 3 ) DM profile in the R3 ROI. The discontinuity in the expected and observed limit in this ROI around 
1 GeV is the result of using only PSF3 type events. See Sec. |III| for more information. The upper right panel is for the Einasto 
profile in the R16 ROI. The lower left panel is the NFW DM profile in the R41 ROI, and finally the lower right panel is the 
Isothermal DM profile in the R90 ROI. Yellow (green) bands show the 68% (95%) expected containments derived from 1000 
no-DM MC simulations (see Sec. |VB|. The black dashed lines show the median expected limits from those simulations. Also 
shown are the limits obtained in our 3.7-year line search [19] and our 5.2-year line search when the assumed DM profiles 
were the same. 


The LAT consists of 16 towers, each includes a tracker module and a calorimeter module |53|. Pass 8 includes 
important updates to the energy reconstruction near the edges of the calorimeter modules (<60 mm from the center of 
the gap) IMlISSj. Events that deposit the majority of their energy (or have their reconstructed centroid) near the edge 
of a calorimeter module are more difficult to reconstruct accurately because of energy leakage of the shower into the 
gaps between modules, or towers. Pass 8 applies an improved handling of this leakage in the energy reconstruction 
algorithms. We show in Fig. [T^the distance of each reconstructed centroid from the center of the calorimeter gap for 
the events passing the comparison selection outlined above. Each calorimeter crystal has a width of 326 mm and the 
gap between modules of 44 mm m- This yields a total width of 370 mm. In this hgure, 0 mm marks the distance 
from the middle of the gap between sets of crystals. The figure at the top also includes a cartoon to illustrate the 
location of the edge of the calorimeter crystal with the center located at 185 mm. 


About half of the overlapping events between Pass 7REP and Pass 8 in the 120-150 GeV energy range were 
reconstructed with centroids near the edges of the towers (<60 mm from the center of the gap). As a consequence, 
these events had the largest differences in reconstructed energy and comprised the tails of the distribution shown on 
the left in Fig. There appears to be a slight enhancement of events where much of the shower was lost between 
modules in the energy range around 133 GeV relative to all events above 20 GeV. 
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FIG. 9. 95% CL tdm lower limits assuming and NFW profile in R180. Yellow (green) bands show the 68% (95%) expected 
containments derived from 1000 no-DM MC simulations (see Sec. VB |. The black dashed lines show the median expected limits 
from those simulations. Also shown are the limits obtained in our 3.7-year line search m and our 5.2-year line search |22| 




(E'(P8)-E'(P7))/E'(P7) Tower Edge to Centroid (mm) 

FIG. 10. Comparison of events selected in R3 from the full data set (5.8 years) passing the Clean event selections simultaneously 
in both Pass 7REP and Pass 8. The left panel shows the fractional energy difference between Pass 7REP and Pass 8. The 
black line shows all events that have a reconstructed energy {E') above 20 GeV. The dashed red line indicates events which 
are with Pass 7REP energies in the range 120-150 GeV. The right panel shows the distribution of the distance from the center 
of the gap of the calorimeter module to the centroid of the cluster. The black curve shows all events above 20 GeV and the 
red curve has the subset in the Pass 7REP range 120-150 GeV. These events are separated by those near the edge of the tower 
(<60 mm from the center of the gap), and those away from the edge of the tower (>60 mm from the center of the gap). The 
blue dashed line marks the edge of the calorimeter module, which is also illustrated at the top of the figure. 


B. Feature in R3 


To understand the impact of Pass 8 on R3, we first considered the same time and energy range as our previous 
3.7-year search m- The feature in Pass 7REP, which was narrower than the energy resolution of the LAT and had a 
local significance of 3.3 (t. With Pass 8 the excess present in Pass 7REP data is reduced to a local significance of 2cr 
as is shown in Fig. El 

We then considered the data for the full 5.8-year time range. Figure [T^ shows the fit for a y-ray line at 133 GeV 
in the 5.8 year Pass 7REP and Pass 8. The Pass 8 data are fit using the method described in Sec. [V| and the Pass 
7REP data are fit similarly, but using the ‘2D’ Deg model described in Sec. IV of Ref. [19]. The Pass 7REP curve 
shows a clear decrease in local significance (from 3.3(T to 2a) with respect to the previous line analysis over a shorter 
time interval [TSj. Similarly for the Pass 8 data, the local significance also decreases {2a to <lcr) using the full 5.8 
year dataset. 
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FIG. 11. Fit at 133 GeV for a 7-ray in the 3.7-year Pass 8 data set using the 2D energy dispersion model in R3. The solid 
curve shows signal and background fitting procedure described in Sec. |V A[ The blue dotted line is the signal model that best 
fits the data. The gray line, which is mostly hidden by the solid curve, is the best fit background. The bin size is such that the 
energy resolution is sampled with 3 bins. 
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FIG. 12. Fit at 133 GeV for a 7-ray in a 5.8-year Pass 7REP (left) and Pass 8 (right) data sets using the 2D energy dispersion 
model in R3. The solid curve shows signal and background fitting procedure described in Sec. |V A| The blue dotted line is the 
signal that best fits the data. The gray line, which is mostly hidden by the solid curve, is the best fit background. The bin size 
is such that the energy resolution is sampled with 3 bins. 


C. Feature in the Earth Limb 


The 7 -ray spectrum of the Earth Limb (see Tab. |T]) is expected to be featureless; however, in the Pass 7REP data 
a 2a feature was found at the same energy as the feature in R3 HUM]. This was a strong indication that the 
feature seen in R3 could have been, in part, a systematic effect. We carried out additional studies with Pass 8 event 
reconstruction and the full dataset to further understand this feature in the Limb. Figure [T^ shows a fit to a 7 -ray line 
at 133 GeV using the full 5.8 year Pass 7REP and Pass 8. We find a slight detection of a line-like feature in both 
Pass 7REP and Pass 8 with a similar fractional size. With Pass 8 the significance increases slightly due mainly to 
the increase in the number of events from the greater acceptance of Pass 8. 


We note that no feature at 133 GeV is present in the GP control region (Sec. III). To try to understand the nature 


of the slight excess in the Earth Limb with no detection in the GP, events in the GP were reweighted in 9 and in 
azimuthal angle, (p, to the distribution in the Limb. This would indicate a dependence of the feature on the particular 
distribution of arrival directions of the 7 rays in instrument coordinates. The reweighting also yielded no detection of 
a line-like feature at 133 GeV in the GP. Additionally the Limb selection criteria was modified (in both 9r and 9^) to 
see if the feature was enhanced or decreased in any particular part of phase space. The only significant change came 
when splitting the Limb data by the signed value of the rocking angle, 9r. The feature appears more significantly 
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FIG. 13. Fit at 133 GeV for a 7-ray in a 5.8-year Pass 7REP (left) and Pass 8 (right) data sets using the 2D energy dispersion 
model in the Earth’s Limb. The solid curve shows signal and background fitting procedure described in Sec. |V A[ The blue 
dotted line is the signal that best fits the data. The gray line, which is mostly hidden by the solid curve, is the best fit 
background. The bin size is such that the energy resolution is sampled with 3 bins. 


(2.6(t) for time intervals when the rocking angle of the LAT was positive {0r >52°), and almost disappears (0.75cr) 
during time intervals when the rocking angle was negative {6^ < —52°). Requiring larger values of \0r\, however, does 
not significantly change the fractional signal or the significance. 


IX. SUMMARY 

In this work, we have presented an updated search for 7 -ray spectral lines using techniques developed in our previous 
line searches [nHnni] across more than three decades of energy using data reprocessed and selected with Pass 8. 
We searched for spectral lines in the energy range from 200 MeV to 500 GeV in five ROIs optimized for signals 
originating from Galactic DM annihilation or decay. We do not find any significant spectral lines and therefore set 
monoenergetic flux upper limits in each of our ROIs. 

Our search improves on our most recent 3.7-year analysis m and 5.2-year analysis |22] with increased exposure and 
a broader energy range. This results in a corresponding increase in the exposure for the entire 5.8-year time range. 
Additionally, the exposure toward our smallest ROIs (R3 and R16) benefited from Fermi operating in a modified 
observing mode that roughly doubled the rate of accumulation of exposure in the GG region for the last five months 
of the dataset relative to normal survey mode. 

Additionally, we have improved our treatment of the systematic uncertainties in this analysis. We explicitly incor¬ 
porate the systematic uncertainties in our limits by including a nuisance parameter in our likelihood function using 
an Sfsyst value based on fits for line-like features along the GP. Our current limits improve, on average, relative to our 
5.2-year limits in part because we use a more realistic Sf^yst value. In fits for which systematic uncertainties were the 
dominant uncertainty (i.e. lower energy fits in larger ROIs that had more than ^10,000 event^, our 5.8-year limits 
are more robust. 

We investigated in particular a previously-reported line-like feature at 133 GeV after an initial excess was found 
with local significance 3.3(T m- Two additional years of data and new event reconstruction and selection algorithms 
see the feature decrease from /(133 GeV)R 3 = 0.61 in the 3.7-year Pass 7REP data set to /(133 GeV)R 3 = 0.07 in 
the 5.8-year Pass 8 dataset. This fractional signal in Pass 8 is consistent with what is also seen in the Earth Limb 
control region. With the entire 5.8-year Pass 8 data set, the local significance has dropped to 0.72a (from 2. Oct in the 
3.7-year data set), which is consistent with most of the original feature originating from a statistical fluctuation. The 
fact that the feature at 133 GeV is still marginally significant in the Earth Limb suggests a small systematic effect at 
this energy; however, no such feature is present in the GP control data set. 

The sensitivity of future line searches with the LAT will increase with continued exposure. Additional improve¬ 
ments require a more sophisticated modeling of the standard astrophysical backgrounds beyond our simple power-law 
approximation. 


For example, in our search, fits with E.y< 10 GeV and Ej< 150 GeV have 5fByBt> <5/atat in R3 and R180 respectively. 
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Appendix A: OTHER CONTROL SAMPLES 

Our main quantification of Sfsyst was determined based on fits in 10° x 10° ROIs along the GP. This is because the 
background 7 -ray emission in the GP is broadly similar to that in our signal ROIs. Systematically-induced line-like 
features in the GP could be produced by approximating the background spectrum as a power law and unmodeled 
energy-dependent variations in the exposure. 

One can also use other control samples to further study systematically-induced line-like features. The two additional 
control samples we consider are the Earth’s Limb and the Vela Pulsar. Both are extensively used by the LAT 
Gollaboration for studies of systematic uncertainties [53] and were also used as control samples in previous line 
searches It should be noted that the intrinsic 7 -ray emission in these two control regions is different from 

that in the GP. 

The Earth’s Limb emission is composed of 7 rays produced via CR interactions in the Earth’s upper atmosphere. 
The Limb dataset is composed of events which pass the selection > 52° & 111° < 9z < 113° (see Tab. |l]). Above 
a few GeV, the energy spectrum of the Limb is well modeled as a power law [3B] . Therefore we expect the power-law 
approximation of background spectra to be better than in the GP. Below a few GeV, modeling the Limb spectrum 
becomes complicated due to the geomagnetic cutoff, so we only fit for line-like features in the Limb with E^> 5 GeV. 

We select a 20° ROI around the Vela pulsar using the ‘Celestial’ data selections outlined in Ta b. In our Vela 
analysis, we use pulsar phases calculated with the Tempo2 packag mm and the standard ephemeri^^ For our line 
fits, we select the on-pulse data (7 rays with phases in the ranges [0^.1,0.3] U [0.5, 0.6]) and model the background in 
energy bin i as a sum of the off-pulse (7 rays with phases in the range [0.7,1.0]) spectrum and an exponential cutoff 
model 


Gbkg.veia = a f — j exp[-{E/Ec)]£{E), (Al) 

where a is the normalization factor, Eq is set to 100 MeV. We fix Ac = 3 GeV (which is slightly different from those 
cited in Ref. [3B]) and let P float free in the fit. Since the Vela spectrum cuts off steeply above a few GeV, we only 
fit for lines up to Ey = 10 GeV, making this region complementary to the higher-energy studies performed using the 
Limb data. The observed fractional signal values in the GP, Vela, and the Limb are shown in Fig. [14] 



FIG. 14. (right) Fractional signals observed in the the Galactic plane (GP) scan (grey dots), from the Vela pulsar (blue squares), 
and from the Earth’s Limb (red triangles). The error bars for the Vela and Limb points are 5/stat = l/V^eff- 

Note that the features observed in the Limb, Vela, and the GP do not have the same energy dependence. This 
is because the inaccuracies of the background model are different in each control region and these errors induced 
different line-like features, especially from the modeling of the background spectra. For example, a slight increase 
of / values is observed in the GP around 1.5 GeV, and a similar increase is observed in the Vela scans, but closer 
to 800 MeV instead. Both background models approximate the energy-dependent exposure similarly, but the flux 


® http://www.atnf.csiro.au/research/pulsar/tenipo2/ 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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models differ since the intrinsic 7 -ray emission in the GP does not cut off like the Vela emission does. Therefore, the 
deficiencies in the total background spectral modeling induce / values with different energy dependencies. However, 
we note that the magnitude of the induced fractional signals is similar across all control regions. This motivated our 
choice of modeling Sfsyst with a general Gaussian envelope centered on zero instead of using a more sophisticated 
energy-dependent model. 


Appendix B: 95% CONFIDENCE LEVEL LIMITS 


We present the 95% CL flux upper limits derived for each of our 5 ROIs in Tabs. [lTHlVl We also give the annihilation 
cross section upper limits for the DM profiles in each ROI where sensitivity to that DM model has been optimized: 
R3 (contracted NFW profile), R16 (Einasto profile), R41 (NEW profile) and R90 (Isothermal profile), and the decay 
lifetime lower limit for R180 (NEW profile). 
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TABLE II. 95% confidence level limits from all ROIs for fit energies from 0.214-5.22 GeV. The first column for each ROI is 
the (‘h-,-,) upper limit in 10“® crn^s”^. The second column for each is the upper limit on {av)-yy in 10“®® cm®s“® for the DM 
profile for which that ROI is optimal. For R180, we state lower limit on r.y^(NFW) in 10®* s. Note that for the energy is 
m^/2. 




R3 


R16 


R41 


R90 


R180 

Energy 

®77 


*^77 

((JU)77 

$77 


$77 


^11 

rVyiy 

(GeV) 


NFWc 


Ein 


NFW 


Iso 


NFW 

0.214 

18.4 

0.766 

45.5 

2.79 

135 

8.52 

315 

25.8 

587 

0.788 

0.234 

15.7 

0.778 

38.0 

2.79 

117 

8.78 

280 

27.4 

530 

0.798 

0.255 

14.8 

0.873 

34.0 

2.97 

102 

9.14 

248 

28.8 

477 

0.813 

0.278 

14.8 

1.03 

32.4 

3.35 

96.4 

10.2 

228 

31.5 

440 

0.810 

0.303 

13.2 

1.10 

30.9 

3.79 

92.0 

11.6 

210 

34.4 

399 

0.820 

0.329 

12.8 

1.26 

28.9 

4.20 

85.5 

12.7 

189 

36.6 

353 

0.852 

0.358 

11.4 

1.32 

26.3 

4.51 

82.4 

14.5 

168 

38.4 

310 

0.893 

0.388 

10.1 

1.37 

21.5 

4.33 

65.5 

13.5 

141 

37.9 

264 

0.965 

0.421 

9.16 

1.47 

18.8 

4.47 

56.5 

13.7 

121 

38.2 

228 

1.03 

0.456 

8.59 

1.61 

16.7 

4.66 

48.9 

13.9 

105 

39.2 

199 

1.09 

0.493 

7.42 

1.63 

14.7 

4.78 

42.2 

14.1 

89.2 

38.8 

171 

1.18 

0.533 

6.45 

1.66 

13.5 

5.13 

37.5 

14.6 

78.8 

40.0 

149 

1.25 

0.576 

5.71 

1.71 

11.8 

5.23 

33.8 

15.3 

70.7 

41.8 

134 

1.28 

0.620 

4.66 

1.62 

10.4 

5.35 

30.6 

16.2 

63.7 

43.8 

121 

1.32 

0.668 

3.83 

1.54 

9.84 

5.88 

29.5 

18.1 

62.0 

49.3 

115 

1.29 

0.718 

3.26 

1.52 

8.78 

6.06 

26.8 

19.0 

57.4 

52.8 

107 

1.29 

0.770 

3.30 

1.77 

7.80 

6.20 

23.1 

18.8 

49.3 

52.3 

76.3 

1.69 

0.826 

2.85 

1.76 

6.74 

6.16 

19.6 

18.3 

39.8 

48.4 

73.4 

1.64 

0.885 

2.26 

1.60 

6.36 

6.66 

17.6 

18.9 

33.6 

47.0 

61.3 

1.83 

0.947 

2.52 

2.05 

7.32 

8.79 

19.3 

23.7 

37.1 

59.5 

68.7 

1.52 

1.01 

2.25 

2.09 

8.68 

11.9 

22.6 

31.8 

44.8 

82.1 

81.9 

1.19 

1.08 

2.45 

2.60 

9.29 

14.6 

24.6 

39.6 

50.5 

106 

92.1 

0.993 

1.16 

2.22 

2.69 

8.68 

15.6 

23.3 

42.8 

49.0 

117 

89.5 

0.956 

1.24 

1.82 

2.51 

7.55 

15.4 

20.8 

43.7 

43.8 

119 

79.3 

1.01 

1.32 

1.47 

2.31 

6.39 

14.9 

18.0 

42.9 

37.0 

115 

67.1 

1.12 

1.41 

1.29 

2.30 

5.12 

13.5 

14.6 

39.6 

29.7 

105 

53.8 

1.31 

1.50 

1.21 

2.45 

4.03 

12.1 

11.5 

35.4 

23.1 

92.3 

41.3 

1.60 

1.60 

0.936 

2.16 

3.17 

10.8 

8.66 

30.3 

17.6 

79.8 

31.9 

1.95 

1.70 

0.815 

2.13 

2.92 

11.3 

7.12 

28.3 

14.2 

73.2 

25.6 

2.27 

1.81 

0.596 

1.77 

2.20 

9.68 

5.66 

25.5 

11.1 

65.1 

20.5 

2.67 

1.93 

0.430 

1.45 

1.84 

9.18 

4.70 

24.1 

9.98 

66.4 

17.9 

2.87 

2.06 

0.428 

1.64 

1.93 

11.0 

4.68 

27.2 

9.70 

73.3 

17.0 

2.84 

2.19 

0.461 

2.00 

1.87 

12.1 

4.63 

30.5 

8.91 

76.4 

15.7 

2.88 

2.33 

0.506 

2.50 

1.75 

12.8 

4.23 

31.6 

7.78 

75.8 

13.7 

3.09 

2.49 

0.365 

2.04 

1.35 

11.2 

3.47 

29.4 

6.88 

76.0 

12.0 

3.32 

2.65 

0.234 

1.49 

1.07 

10.1 

2.87 

27.6 

6.19 

77.5 

10.7 

3.51 

2.82 

0.239 

1.72 

0.844 

8.99 

2.48 

27.1 

5.61 

79.8 

9.70 

3.62 

3.00 

0.282 

2.31 

0.738 

8.92 

2.31 

28.6 

5.26 

84.8 

8.63 

3.82 

3.20 

0.222 

2.05 

0.780 

10.7 

2.22 

31.1 

4.28 

78.1 

7.30 

4.25 

3.40 

0.289 

3.03 

0.915 

14.2 

2.07 

32.9 

3.85 

79.5 

6.47 

4.50 

3.62 

0.326 

3.86 

1.02 

17.9 

2.11 

38.0 

3.22 

75.3 

5.76 

4.75 

3.85 

0.321 

4.31 

0.925 

18.4 

1.95 

39.7 

2.83 

74.9 

4.91 

5.24 

4.09 

0.271 

4.11 

0.764 

17.2 

1.52 

35.0 

2.43 

72.7 

4.01 

6.03 

4.35 

0.175 

3.00 

0.715 

18.1 

1.25 

32.5 

2.22 

75.1 

3.66 

6.21 

4.63 

0.180 

3.48 

0.561 

16.0 

1.11 

32.6 

2.00 

76.4 

3.15 

6.80 

4.91 

0.157 

3.43 

0.445 

14.4 

0.964 

31.9 

1.68 

72.5 

2.84 

7.10 

5.22 

0.115 

2.82 

0.263 

9.59 

0.656 

24.5 

1.33 

64.7 

2.35 

8.10 
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TABLE III. 95% confidence level limits from all ROIs for fit energies from 5.54-43.8 GeV. The first column for each ROI is 
the ('h.y.y) upper limit in 10“^° cm^s“^. The second column for each is the upper limit on {av)-yy in 10“^® cm®s“^ for the DM 
profile for which that ROI is optimal. For R180, we state lower limit on r.y^(NFW) in 10®* s. Note that for the energy is 
m^/2. 




R3 


R16 


R41 


R90 


R180 

Energy 

^77 


*^77 

{av)jj 

"^77 


"^77 

((JU)77 

"^77 

TVyiy 

(GeV) 


NFWc 


Ein 


NEW 


Iso 


NEW 

5.54 

4.58 

0.127 

21.6 

0.888 

64.1 

2.69 

no 

6.04 

196 

9.13 

5.87 

4.62 

0.144 

20.6 

0.950 

63.7 

3.01 

102 

6.27 

187 

9.03 

6.23 

5.65 

0.198 

17.0 

0.881 

61.0 

3.25 

115 

7.96 

220 

7.23 

6.60 

7.84 

0.309 

22.5 

1.31 

75.0 

4.48 

144 

11.2 

263 

5.70 

6.99 

4.74 

0.209 

26.0 

1.70 

74.2 

4.98 

160 

13.9 

271 

5.23 

7.40 

5.43 

0.269 

33.4 

2.45 

76.2 

5.73 

152 

14.8 

259 

5.16 

7.83 

9.21 

0.510 

25.3 

2.08 

55.0 

4.62 

119 

13.1 

211 

6.01 

8.28 

6.14 

0.381 

18.6 

1.71 

39.8 

3.74 

100 

12.3 

177 

6.76 

8.76 

4.63 

0.321 

22.1 

2.27 

37.5 

3.95 

91.0 

12.5 

153 

7.37 

9.26 

3.62 

0.281 

12.7 

1.46 

28.8 

3.39 

74.2 

11.4 

131 

8.19 

9.79 

3.42 

0.297 

7.59 

0.974 

33.0 

4.35 

61.7 

10.6 

115 

8.81 

10.4 

5.61 

0.543 

8.14 

1.17 

28.9 

4.25 

59.6 

11.4 

106 

8.99 

10.9 

6.51 

0.705 

12.4 

1.98 

35.3 

5.80 

78.9 

16.9 

121 

7.48 

11.6 

6.41 

0.775 

16.0 

2.87 

36.4 

6.67 

73.9 

17.6 

117 

7.32 

12.2 

7.26 

0.980 

7.72 

1.54 

21.3 

4.35 

59.9 

16.0 

111 

7.28 

12.9 

5.15 

0.775 

7.83 

1.75 

23.2 

5.30 

56.7 

16.9 

111 

6.91 

13.6 

2.97 

0.499 

8.42 

2.09 

21.9 

5.57 

41.3 

13.7 

84.1 

8.64 

14.4 

2.86 

0.536 

5.91 

1.64 

15.3 

4.34 

27.4 

10.2 

60.9 

11.3 

15.2 

3.00 

0.626 

6.52 

2.02 

15.8 

5.01 

30.1 

12.4 

61.4 

10.6 

16.1 

2.76 

0.645 

6.66 

2.30 

17.1 

6.06 

40.6 

18.7 

68.4 

9.02 

17.0 

4.07 

1.06 

8.18 

3.15 

23.0 

9.08 

39.8 

20.4 

68.4 

8.54 

17.9 

2.50 

0.726 

10.2 

4.37 

20.8 

9.19 

42.5 

24.4 

62.6 

8.83 

18.9 

3.40 

1.10 

5.79 

2.78 

13.8 

6.79 

34.9 

22.4 

50.2 

10.4 

20.0 

5.54 

2.00 

3.65 

1.96 

11.3 

6.21 

19.6 

14.0 

32.5 

15.3 

21.1 

4.14 

1.67 

6.56 

3.92 

11.2 

6.84 

12.2 

9.74 

22.3 

21.0 

22.3 

2.10 

0.946 

3.66 

2.44 

5.57 

3.81 

9.01 

8.02 

18.7 

23.7 

23.6 

1.90 

0.957 

3.74 

2.78 

6.43 

4.91 

9.89 

9.83 

17.7 

23.7 

24.9 

1.47 

0.829 

2.97 

2.47 

8.70 

7.43 

14.8 

16.4 

24.2 

16.4 

26.4 

1.41 

0.888 

3.78 

3.51 

11.4 

10.9 

19.7 

24.5 

38.1 

9.86 

27.9 

1.77 

1.24 

4.56 

4.74 

14.4 

15.3 

29.8 

41.3 

46.3 

7.67 

29.5 

1.22 

0.958 

7.05 

8.19 

15.7 

18.7 

26.4 

40.9 

44.9 

7.50 

31.2 

0.753 

0.661 

4.37 

5.68 

9.49 

12.6 

20.5 

35.5 

37.5 

8.47 

33.0 

0.708 

0.695 

3.28 

4.77 

4.12 

6.14 

13.7 

26.6 

31.2 

9.62 

34.9 

1.29 

1.42 

4.17 

6.79 

5.95 

9.93 

14.9 

32.4 

25.4 

11.2 

36.9 

2.51 

3.08 

4.71 

8.58 

10.6 

19.9 

23.2 

56.4 

28.8 

9.34 

39.0 

2.81 

3.87 

3.18 

6.48 

9.47 

19.8 

18.4 

50.1 

25.4 

10.0 

41.3 

2.45 

3.79 

3.07 

7.01 

7.00 

16.4 

14.0 

42.8 

23.1 

10.4 

43.8 

3.70 

6.41 

4.71 

12.1 

7.13 

18.8 

12.7 

43.6 

20.7 

11.0 
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TABLE IV. 95% confidence level limits from all ROIs for fit energies from 46.4-462 GeV. The first column for each ROI is the 
(i&.y.y) upper limit in 10“^® crn^s”^. The second column for each is the upper limit on {av)jry in 10“^® cm®s“^ for the DM 
profile for which that ROI is optimal. For R180, we state lower limit on r.y^(NFW) in 10®* s. Note that for the energy is 
m^/2. 




R3 


R16 


R41 


R90 


R180 

Energy 


{ov)^^ 

^77 


"^77 

((JU)77 

$77 


$77 

TTyiy 

(GeV) 


NFWc 


Ein 


NEW 


Iso 


NEW 

46.4 

3.37 

6.54 

5.66 

16.3 

7.64 

22.5 

11.6 

44.7 

16.4 

13.0 

49.1 

2.58 

5.62 

6.40 

20.7 

5.66 

18.7 

6.18 

26.6 

10.9 

18.4 

52.1 

1.27 

3.12 

4.56 

16.6 

3.07 

11.4 

4.55 

22.0 

9.86 

19.3 

55.2 

1.23 

3.38 

3.96 

16.2 

3.88 

16.2 

4.73 

25.7 

10.2 

17.7 

58.6 

2.30 

7.13 

4.85 

22.3 

5.87 

27.7 

7.28 

44.6 

10.0 

16.9 

62.2 

1.99 

6.95 

3.32 

17.2 

2.85 

15.1 

4.38 

30.2 

5.77 

27.6 

66.0 

1.17 

4.59 

1.82 

10.6 

2.11 

12.6 

5.03 

39.1 

5.38 

27.9 

70.1 

1.17 

5.18 

1.90 

12.5 

4.49 

30.2 

5.92 

51.9 

5.61 

25.2 

74.5 

1.11 

5.56 

3.63 

26.9 

6.61 

50.3 

4.81 

47.6 

9.46 

14.1 

79.2 

0.543 

3.08 

1.48 

12.4 

3.31 

28.5 

4.27 

47.8 

8.32 

15.0 

84.2 

0.448 

2.87 

0.951 

9.03 

1.88 

18.3 

2.82 

35.7 

4.95 

23.8 

89.6 

0.396 

2.87 

0.947 

10.2 

2.36 

26.0 

4.20 

60.2 

6.77 

16.3 

95.4 

0.343 

2.82 

0.891 

10.8 

2.88 

35.9 

7.63 

124 

9.86 

10.5 

102 

0.619 

5.77 

2.29 

31.6 

4.25 

60.0 

6.57 

121 

9.60 

10.2 

108 

0.541 

5.73 

4.89 

76.5 
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